Berry Phase in Atom-Molecule Conversion Systems and Fractional Monopole 
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We investigate the geometric pliase or Berry phase of adiabatic quantum evolution in an atom- 
molecule conversion system, and find that the Berry phase in such system consists of two parts: 
the usual Berry connection term and a novel term from the nonlinearity brought forth by the 
atom-molecule conversion. The geometric phase can be viewed as the flux of the magnetic field 
of a monopole through the surface enclosed by a closed path in parameter space. The charge of 
the monopole, however, is found to be one third of the elementary charge of the usual quantized 
monopole. 
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Berry phase [H , which reveals the gauge structure asso- 
ciated with a phase shift in adiabatic processes in quan- 
tum mechanics, has attracted great interest in physics Ql- 
One classical example of Berry phase is a spin half par- 
ticle aligned to a magnetic field, and the field is made to 
rotate adiabatically in a 3-D parameter space (see Fig.l 
a) . The Berry phase of such system has been interpreted 
as the flux of a magnetic field of a quantized monopole 
through the surface enclosed by the loop in parameter 
space. 

On the other hand, association of ultracold atoms 
into molecules is currently a topic of much experimen- 
tal and theoretical interest Q with important applica- 
tions ranging from the search for the permanent electric 
dipole moment [3| to BCS-BEC (Bose-Einstein conden- 
sate) crossover physics [Sl . Through Feshbach resonance 
or photoassociation [3| , a pair of atoms can convert 
into a bounded molecule. The atom-molecule conversion 
under mean field treatment is governed by a nonlinear 
Schrodinger equation, in which the nonlinearity is from 
the fact that two or more atoms are needed to form one 
moleculeQ. Since the adiabatic manipulation is an opti- 
mal way to yield high conversion efficiency, great efforts 
and big progress have been made towards the adiabatic 
condition and adiabaticitv [lo| of the nonlinear quan- 
tum evolution. Nevertheless, the knowledge of the Berry 
phase for the adiabatic evolution in such nonlinear sys- 
tem is very limited. This system not only lacks super- 
position principle due to the presence of nonlinearity [U 
but also have no ?7(l)-invariance because the chemical 
potentials of atomic component and molecular compo- 
nent are not identical |il2i]. 

In this letter we formulate the adiabatic geometric 
phase in a general formulism for the atom-molecule con- 
version systems and derive the explicit expression of the 
Berry phase analytically. We find strikingly that the cir- 
cuit integral of Berry connection of the instantaneous 
eigenstate alone can not account for the geometric phase, 
while a novel term due to the nonlinearity brought forth 
by the atom-molecule coupling emerges. Only with the 
inclusion of this additional contribution, the total geo- 



metric phase can be interpreted as a flux of a magnetic 
fleld of a monopole through the surface enclosed by the 
closed path in parameter space. There exists a forbidden 
cone in Berry sphere (i.e., parameter space, see Fig. 1 
a') and the Bloch sphere representing projective Hilbert 
space is deformed dramatically (Fig.l b'). The charge of 
the monopole is found to be one third of the elementary 
charge of the usual quantized monopole. 

Let us consider an atom-molecule system with energy 
7i(V', V"^; R);!!! which R denotes all the system parame- 
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FIG. 1: (color online) The parameter space (a) and Bloch 
sphere (b) for a spin half particle in magnetic field, (a') 
and (b') are the parameter space and Bloch sphere for the 
atom-molecule conversion system, respectively. The param- 
eters change adiabatically along a close path shown as the 
green circles in the parameter spaces or Berry spheres . Ac- 
cordingly, the eigenstate will evolve and form a close path 
schematically plotted as green circles on the Bloch spheres. 
The gray cone in (a') is the boundary for which 9 — 2tt/3., 
inside which, i.e., 9 > 27r/3, no eigenstate exists. See text for 
details. 
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ters that vary in time slowly, 'ip = {V'i} and 'ip'^ = {V']} 
(with i = 1, • ■ ■ ,N) are the annihilation and creation 
field operators for atom with i — I,-- - ,M and for 
molecule with i = M + 1, • • • ,N. They obey the com- 



mutation relations 



Sij for bosons. 



Under 



the mean field treatment, tp and are replaced by 
complex numbers ^p and ip*. It is convenient to write 
n{ip,'p*;R) = J2t,j^tTzj{^,^*;'R)i'j, then with the 
help of the above commutation relations, we obtain fol- 
lowing nonlinear Schrodinger equations (?i = 1) 



k 



dt 



in which the Hamiltonian 



V^* ; R) = T,fe (V, V* ; R) + ^ V' ■ 



or. 



(1) 



(2) 



The above system is not invariant under usual U{1) 
transformation, instead, it is invariant under the follow- 
ing co-diagonal U{1) transformation. 



U{r]) = e*®("\e(??) 



Vl'' , 
K7?/^-^^ )' ^'^> 



in which is the K rank unit matrix and k is an in- 
teger related with molecule structure, e.g. for diatomic 
molecule k = 2 1^. Obviously, for k = 1 the system re- 



duces to an atomic system and U{ri) is just the ordinary 
U{1) transformation. 

The eigenequation of the above system is 

J2 Hjk (0(R) , ^* (R) ; R) 0, (R) = A*, (R)^, (R) , (4) 



where /ij(R) = m(R) for j < M (atom) and /ij(R) = 
K/i(R) for j > M(molecule). 

The above eigenequation defines the eigenfunction (p 
and the eigenvalue (or chemical potential) n that are 
functions of the adiabatic parameter R. 

The equation ^ and its conjugate construct a canon- 
ical structure of classical dynamics with the energy 
H{ip,'p*;'R) as classical Hamiltonian and {ip,iip*) as a 
canonical variable pair. The gauge symmetry of Ti given 
by ([3]) implies that the total atoms number is conserved 
and the dynamics of the overall phase can be separated 
from the rest of the degrees of freedom [l^. For simplic- 
ity and without losing generality we denoteA = argV'i 
and set total phase as A for atomic components and 
kX for molecular components, respectively. The normal- 
ization condition is Y^fLi IV'iP + i^Y.^=m+i I^jI^ = 1- 
The other variables form a close set of Hamiltonian 
dynamics. We choose a new set of canonical vari- 
ables {q,p) as q = {qi,- ■ ■ ,qi,- ■ ■ ,qN~i) with qt = 
- arg(V'i+i) + Ki+i arg(-0i), p = (pi, ■■■ ,Pi,--- ,Pn-i) 
with Pi — l^j+ip, in which Ki — 1 for i < M and Ki = k 



for i > M. From ([T]) and its canonical structure and us- 
ing the normalization condition, we obtain the dynamical 
equations for the overall phase and other variables. 



d\ dq 
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Hp, q), 



p = 

in which A(p, q) 

gi(arg(i/jfc)-arg(-0i)) jv^ 
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Re 



m 

dp ' 



(5) 
(6) 



Pk] and Tikip,q) 



and Tij are the 



For a linear quantum case, both Hij 
functions of the parameter R only, so that the last term 
in Eq.(5) vanishes, i.e., A{p,q) — 0. The second term in 
the right-hand of Eq.(5) is the energy, whose time inte- 
gral gives so called dynamical phase. The time integral 
of the first term is the Aharonov-Anandan phase for a 
cyclic quantum evolution 1J| . The above observation is 
readily extended to the adiabatic evolution of a quantum 
eigenstate, because the adiabatic theorem of quantum 
mechanics dictates that an initial nondegenerate eigen- 
state remains to be an instantaneous eigenstate and the 
evolution will be cyclic when the parameters move slowly 
in a circuit. In this case, the second term is the eigenen- 
ergy and the first term is just the Berry connection, i.e. 
i < ^(R)|V|0(R) >. Then the Berry phase equals to the 
circuit integral of the Berry connection. 

However, for our nonlinear system, the contribution of 
the last term in (O should be taken into account. Notice 
that the chemical potential is usually not identical to 
the energy while the dynamic phase should be the time 
integral of the chemical potential, we need to evaluate 
following quantity in adiabatic limit. 



S(p, q; R) - n{p, q) + A{p, q) - m(R), 



(7) 



We denote p = p(R) + dp and q = g(R) -I- Sq . Here p(R) 
and g(R) are the fixed points of Eq.(6) that are local en- 
ergy minima of system and therefore correspond to the 
eigenstates defined by (4) . The vector {Sp, 5q) represents 
the correction to the adiabatic eigenstates in the order of 
I ^ I . As will be shown, {Sp, 6q) contains some secular 
terms in addition to the rapid oscillations. These secu- 
lar terms will be accumulated in the nonlinear adiabatic 
evolution and contribute to the geometric phase. 

We expand the quantity S(p, g;R) around the fixed 
point. 



^{p, q; R) 



OA. 
dp 



\iP,q) 



Sp - 



OA 
dq 



(^p^g^Sq + o{5q\Sp^). (8) 



Here we use the relations: H(p, q) + A{p, q) — ^(R) and 
<9H(p, q)/dp\p,q = dnip, q)/dq\p.q = 

On the other hand, the {Sp, Sq) can be evaluated from 



3 



following Hamiltonian equations, 



P 



dqdp 



Omitting the higher order terms, keeping the secular 
terms by average over the fast oscillations, we obtain 

i{Sp),{Sq)f-n-\^il,^Rr, (11) 



where the matrix 



dqdp 



dqdq 



d^H{p,q) a'n(p,q) 



\ dpdp dpdq / 

the Hessian matrix of the classical Hamiltonian. 

Combining Eq.(8) and (11) and with the help of Eq.(5), 
we find that, except for the dynamical phase (i.e., the 
time integral of the chemical potential), the eigenstate of 
the nonlinear atom-molecule conversion system acquires 
following additional phase during the adiabatic cyclic 
evolution 1 1611 ■ 



la 



dK dA 



n- 



. dp dq J, 



dR.(12) 



In contrast to previous works |17f |, the adiabatic ge- 
ometric phase in the atom-molecule system is dra- 
matically modified. The first term is the usual 
expression of the Berry phase that can be rewrit- 
ten as the circuit integral of the Berry connection 
^ ((0(R)| V0(R)) - (V^(R)| 0(R))). The novel sec- 
ond term indicates that, the high-order correction to adi- 
abatic approximate solution that is negligible in linear 
case, could be accumulated in the nonlinear adiabatic 
evolution with an infinite time duration in adiabatic limit 
and contributes a finite phase with geometric nature. 

As an illustration, we consider following atom-molecule 
conversion model whose energy takes the form. 



n = 



Rcos0 



3 i?sine' 



(^e-"f'iJ;UjUj2 + h.c}j , (13) 



where ip — (V'i,'!/'2) and ■(/'^ = annihila- 
tion and creation operators for atom and molecule respec- 
tively, the terms ■il;\ijj\ip2 + h.c. describe coupling between 
atom pairs and diatomic molecules, and R — {R, 9, ip) are 
parameters. Obviously, the system is invariant under the 

transformation U{i]) ~ e*®(''\0(77) ^ ^ g 2?^ 

Let us rewrite H(V',V'*;R) = J2^J i^tT,j{i^,i;*;R)i;j , 
where the matrix elements Tn = -T22 = T12 = 



3 Rsind, 



'21 Y 8 2 

equation takes the form of Eq.(l) with 



ipl and the nonlinear Schrodingcr 



H(V',V*;R) 



|e-"^i?sin6'V^* 

FI.COS0 



(14) 

where and ■02 are complex amplitudes for 

atom and molecule respectively. The projective 
Hilbert space is spanned by the vector ria = 

(2V2 Re[iri?^2l^V2Imi{ri?i;2),\^i\^ " ^\H^)- Ob- 
viously, every point in this space corresponds to a class 
of quantum states among which the states are only dif- 
ferent in co-diagonal total phases (see Eq. ([3])). With 
the normalization condition iV'iP + 2|'(/'2p = 1; we plot 
the projection space in Fig. 2(b'), which is a "tear-drop' 
shaped surface[18|. 

The eigenequations take the forms of 

tions are solved and the eigenfunctions are obtained as 
follows. 



H 



;R) 



The eigenequa- 



4>2 



-cos6l±l) -± 

01 



iip/2 



1 



2|02^P, 



(15) 



\/6sin6' 

with the eigenvalue (or chemical potential) ^± = 
I (cos6l± 1). 

Following our general formalism, we choose the to- 
tal phase as A = arg-^i and define the canonical pair 
a,s q — — argi/)2 + 2arg0i and p = |V'2p- Hence, we 
have p — 102 P q ~ tp. Substituting the concrete 
expressions of matrix {Tij} to the definition of the quan- 
tity A, we get A(p, q) = ■^§-^^2^(1 - 2p)^cos {q - ip). 
On the other hand, from the system energy 7Y(p, q) — 
^^f^ (1 - 3p) + sin 61(1 - 2p)^cos{q - ^), we get 

1 



-1 _ Vs 



(l+6p) 





After lengthy 



\ (l-2p)VP 

calculation, we obtain the Berry phase according to for- 
mula (fl! 



7g 



pdip 



(1 - 6p)p 
l + 6n 



dip 



^ (j) (l=Fcos6') dip, 



(16) 
(17) 



The above theoretical formulation on the Berry phase 
has been verified numerically by directly integrating the 
Schrodinger equation. On the other hand, we recognize 
that the above system although admits the quantal equa- 
tions of motion, appears formally to have classical struc- 
ture if we regard the total phase and total particle num- 
ber as a pair of canonical conjugate variables. We thus 
could exploit this particular feature to construct a canon- 
ical transformation to action- angle variables. With in- 
cluding the canonical motions in projective Hilbert space 
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represented by Eq. (6) , we have derived the Hannay's an- 
gle of a geometric nature associated with the adiabatic 
evolution. The Hannay's angle is found to exactly equal 
to minus Berry phase of Eq. Td] . This fact indicates a 
novel connection between Berry phase and Hannay's an- 
gle in contrast to the usual derivative formQ, and sup- 
ports our new expression of Berry phase from the other 
aspect. 

For the linear systems, such as the spin-half system, 
i.e., H — — -iR-cr where a are pauli matrices and 
R = (i? sin 6 cos (j), R sin 6 sin 0, R cos 9) is a. vector in the 
3-D parameter space. The Berry phase equals to the cir- 
cuit integral of the Berry connection and is interpreted 
as the flux of the magnetic field of a virtual quantized 
monopole through the surface enclosed by the loop in 
parameter space, i.e., Bm = go-^ with the elementary 
charge 5o = ^- In general, the degeneracies of the spec- 
trum in parameter space play an important role in con- 
nexion with the geometric phase. Each degeneracy can 
be seen as a charge distribution located at the contact 
point between energy surfaces. Because the eigenstates 
are smooth and single valued outside the degeneracies, 
the total charge of the distribution, i.e., the monopole 
charge is necessarily an integer multiple of the elementary 
charge go = 1/2. The non-elementary monopoles with 
integer multiples of go have been found in case of light 



propagating and in condensed matter physics |21l . [2 



The mechanism for the production of monopole charges 
larger than the elementary go is due to constraints that 
act on the system ■ 

For our nonlinear system, when the parameters R = 
{R, 9, (fi) are considered as spherical coordinates of a vec- 
tor in a 3-D space, and then from P7)l we get the vec- 
tor potential, A = ^ ^^RsinP ^v- Here, for convenience 
we only consider the branch of Hence, the Berry 
phase of atom-molecule conversion system can also be 
interpreted as the fiux of a magnetic field of a virtual 
monopole through the surface enclosed by the closed 
path in parameter space (see in fig. 2(a')), that is, 
B = \/ X A.spin — gi^- Strikingly, the monopole charge 
9 — Iso, one third of the elementary charge. We at- 
tribute the fractional charge to the symmetry break- 
ing of the parameter space by the boundary. From 
(fT5|) . we see \ip2\ increases with 9 monotonously. When 
9 — 27r/3, IV'21 — l/v^ reaching its extreme value (since 
iV'iP + 2|'(/'2p = !)■ It implies that there is no eigenstate 
in the regime 9 > 27r/3, i.e., the Berry sphere of this sys- 
tem is a defected sphere with a forbidden cone bounded 
by 9 = 2tt/3. This curious structure has been illustrated 
in fig. l(a'). 

In summary, we have investigated the adiabatic geo- 
metric phase in the atom-molecule conversion systems. 
A novel formula of geometric phase is derived and an ex- 



otic monopole with fractional elementary charge is found. 
The above phase and monopole are expected to be ob- 
served in the future ultracold atom experiments. 
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